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ABSTRACT 

VDB0-B195D is a massive, blue star cluster in M31. It was observed as part of the Beijing- Arizona- 
Taiwan-Connecticut (BATC) Multicolor Sky Survey using 15 intermediate-band filters covering a 
wavelength range of 3000-10,000 A. Based on aperture photometry, we obtain its spectral-energy dis- 
tribution (SED) as defined by the 15 BATC filters. We apply previously established relations between 
the BATC intermediate-band and the Johnson-Cousins UBVRI broad-band systems to convert our 
BATC photometry to the standard system. A detailed comparison shows that our newly derived VRI 
magnitudes are fully consistent with previous results, while our new B magnitude agrees to within 2a. 
In addition, we determine the cluster's age and mass by comparing its SED (from 3000 to 20,000A, 
comprising photometric data in the 15 BATC intermediate bands, optical broad-band BVRI, and 
2MASS near-infrared JHKg data) with theoretical stellar population synthesis models, resulting in 
age and mass determinations of 60.0 ± 8.0 Myr and (1.1 — 1.6) x lO^M©, respectively. This age and 
mass confirms previous suggestions that VDB0-B195D is a young massive cluster in M31. 
Subject headings: galaxies: individual (M31) - galaxies: star clusters - galaxies: stellar content 



1. INTRODUCTION 

Young massive star clusters (YMCs) are among the 
main objects resulting from violent star-forming episodes 
triggered by galaxy collisions, mergers, and close en- 
counters (see lde Grijs &: Parmentieil2007l and references 
therein). They are also referre d to as 'youn g populous 
clusters,' a term first coined bv iHodg^ (|1961| ). who used 
it to describe 23 clusters containin g bright , blue stars 
in the Large Magellanic Cloud. In iHodg^ ()196lD . the 
'young' aspect is demonstrated by the fact that all clus- 
ters have main sequences that extend to absolute magni- 
tudes brighter than My = 0, while 'populous' describes 
their richness (stellar membership). H owever, YMCs are 
also observed in quiescent galaxies (iLarsen fc Richtleil 
11999) and in the disks of isolated spirals, although higher 
cluster-formation efficiencies are associated with environ- 
ments exlnhitnm high star-formation rates (see iLarsen 
120041 iCao fc Wd I2007L and references therein). It has 
become clear that, in many ways, YMCs resemble young 
versions of the old globula r clusters (GCs) a ssociated 
with all large galaxies (see ILarsen "eFaL|[20Ql and ref- 
erences therein). YMCs are seemingly absent in the 
Milky Way; possibly the best example of a Galactic 
YMC is Westerlund 1, a hea vily reddened cluster with 
an age and mass of 4-5 Myr (ICrowther et"alll2006h and 
Mci - 10^ Mq (jClark et al.l '2005\ respectively. 

Since the pioneering work of Tinsley (1968.. 11971 and 
iSearld (|1973| ). evolutionary population synthesis model- 
ing has become a powerful tool to interpret integrated 
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spectrophotometric observations of galax ies and their 
comp onents, such as star clusters (e.g., lAnders et al.l 
l2004f ). The evolution of star clusters is usually modeled 
by means of the simple stellar population (SSP) approx- 
imation. An SSP is defined as a single generation of 
coeval stars formed from the same progenitor molecular 
cloud (thus implying a single metallicity) , and governed 
by a given stellar initial mass function (IMF). 

Age and metallicity are two basic star cluster parame- 
ters. The most direct method to determine a cluster's 
age is by employing main-sequence photometry, since 
the absolute magnitude of the main- sequence turnoff is 
predominantly affected by age (see iPuzia et al.l I2002L 
and references the rein). However, until recently (cf. 
iPerina et al.l l2009f ). this method was only applied to 
the s tar clusters in the Milky Wa y and its satel lites 
(e.g.. lRich eraIll200lD . although Br own et al.l H2004D es- 
timated the age of an M31 GC using extremely deep im- 
ages observed with the Hubble Space Telescope (JfSTj's 
Advanced Camera for Surveys. Generally, the ages of 
extragalactic star clusters are determined by compar- 
ing their observed spectral-energy distributions (SEDs) 
and / or spectroscopy with the p redictions of SSP mod- 
els (IWiUia ms & Hodge 2001a.b: de Griis et al.li200 3a.b/ 



Bik et al.l 1 2003: Jiang et al.l 1200 3: Bcasl ev et al, 
Puzia et al.| .2005; Ma et al.l 12006; .Fan et al 



Ma et al.l 1200^12009 : Caldwell et al.| 120091 : IWang eTaLl 
20101) . Nevertheless, SSP models assume that cluster 



IMFs are fully populated, i.e., that clusters contain in- 
finite numbers of stars with a continuous distribution of 
stellar masses, and that all evolutionary stages are well 
sampled. Real clusters, however, contain a finite number 
of stars. Therefore, a disagreement between the observed 
cluster colors and theoretical col ors derived from SSP 
models may bec ome apparent (see IPiskunov et 311120091 : 
IPopescu fc Hanson 2010. and references therein). Other 
limitations inherent to SSP models arise from our poor 
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understanding of some advanced stellar evolutionary 
stages, such as the supergi ant and the asymptotic-g iant- 
branch (AGB) phases (see lBruzual fc Charloti[2003l and 
references therein). 

Located at a distance of 785 ± 25 kpc, correspond- 
ing to a distance modulus of (m — M)o = 24.47 ± 
0.07 mag (jMcConnachie et al.ll2005[) . M31 is the near- 
est and largest spiral galaxy in the Local Group 
of galaxies. It has been the subject of many GC 
studies and surveys, dating back to the early study 
of Hubble (19 32). Based o n previous publications 
gubbld 119321: [Sevfert fcl^assa u 1945; HiltnCTL^^" " 



iMavall fc EggenI 119531: iKron fc Mavall .1960^ )~ IVetesnikl 
( 19621) compiled the first large M31 GC catalog, con- 
taining UBV photometric data of approximately 300 
GC candidates. Over the past decades, several ma- 
jor catalogs of M31 GCs and GC candidates have 
been p ublished, includi ng major efforts by the B ologna 
group (iBattistini et al.lll 980. 1987, 1993).'Barmb v et al 
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Galleti et al.Tl2004,2005^,200&, 2007,) ..Kim et al 



Caldwell et al.l 1120091 ) ■ and TPeacock et al. 1 (120101) . 



Following on from the first extensive spectroscopic survey 
of M31 GCs by van den Bergh ( 1969) , a signifi cant num- 
ber o f authors (e.g.. iHu chra et al. 1982: H uchra et al.l 
1991| iDubath fc Gr ilhnair 1997; ,F edcrici e t al.l 119931: 



Jablonka et al.l 11998 : Bar mbv et al l 12 000: P errett et al l 
200a IGalleti et al.ll2006t ILee et al.1 12008. and references 



therein) have studied their spatial, kinematic, and chem- 
ical (metallicity) properties. 

M31 is known to host a large number of young 
star cl usters (e.g.. iFusi Pecci et al.l 120051 : iCaldwell et al.l 
20091: IWang et al.l 120101 and references therein). 



Fusi Pecci et al.l ( 20051) presented a comprehensive study 



of 67 very blue star clusters, which they referred to as 
'blue luminous compact clusters' (BLCCs). Since they 
are quite bright (—6.5 < My < —10.0 mag) and very 
young (< 2 Gy r), BLCCs rriay be eguivalent to YMCs 
(see for details jPerina et al.l 200t^ [2 0101. To ascertain 
their properties. iPerina et al. ( 20091 . 1201Q ) performed an 
imaging survey of 20 BLCCs in the disk of M31 using the 
ffST's Wide Field and Planetary Camera-2 (WFPC2). 
They obtained the reddening values, ages, and metallic- 
ities of their sample clusters by comparing the observed 
color-magnitude diagrams (CMDs) and luminosity func- 
tions with theoretical models. 

VDB 0-B195D was first detected by Ivan den Bergh I 
(119691). Its color is extreme ly blue (e.g., U - B = -0.48 
mag: Ivan den Bergh II1969I ) and it is very bright in blu e 
bands (e.g., U = 14.66 mag: Ivan den Bergh I IT969I ). 
As a consequence, Ivan den Bergh I (|1969D asserted that 
VDB0-B195D is the brightest open cluster in M31. He 
determined an integrated stellar spectral type equiva- 
lent to AO, which implies that the cluster contains mas- 
sive stars. In addition, VDB0-B195D is particularly ex- 
tended and most previous photometric studies did not 
include the full extent of the obje ct's light distribution 
(see for details [Perina et al.l 120091 ). We will provide an 
overview of previous studies that included the cluster 
in §2.1. It was observed as part of the galaxy calibra- 
tion program of the Beijing- Arizona- Tai wan-Connecticut 
(BATC) Muhicol or Sky Survey (e.g., iFan et al.l IT996I: 
iZheng et al.lll999l ) in 15 intermediate-band filters. Com- 
bined with photometry in optical broad-band BVRI and 
near-infrared JHKs filters from the Two Micron All Sky 



Survey (2MASS) taken from IPerina et al.l ()2009D . we ob- 
tained the SED of VDB0-B195D in 22 filters, covering 
the wavelengh range from 3000 to 20,000 A. 

In this paper, we describe the details of the observa- 
tions and our approach to the data reduction in §2. In 
§3, we determine the age and mass of VDB0-B195D by 
comparing observational SEDs with population synthesis 
models. We discuss the implications of our results and 
provide a summary in §4. 

2. OPTICAL AND NEAR-INFRARED OBSERVATIONS OF 
THE YMC VDB0-B195D 

2.1. Historical overview 

VDB0-B195D was first given the designation '0' (i.e., 
V DBO), the brightest op e n cluster in M31, in the catalog 
of Ivan den Bergh I (|1969t ). IBattistini et~aLl (|1987[ ) identi- 
fied VDB0-B195p inde pendently and called it B195D. In 
IBattistini et al.l (|1987| ). B195D was given a low level of 
confidence (class D) of being a genuine cluster (classes 
A and B were assigned very high and high levels of 
confidence, respectively). It was only re cently indepen- 
dently confirmed to be a single object. ICaldwell et al.l 
(|2009| ) presented a new catalog containing 670 likely star 
clusters, stars, possible stars, and galaxies in the field 
of M31, all with updated high-quality coordinates ac- 
curate to 0.2", based on images from either th e Local 
Group Galaxies Survey (LGGS) ()Massevl l2006l) or the 
Digitized Sky Survey (DSS). They use the designation 
VDB0-B195D, associated with ao = 00'^40"^29^43 and 
So = -h40°36'14".8 (J2000.0), which are the coordinates 
we ad opt in this paper. Independently, IPerina et aX] 
(|2009D studied the properties of VDB0-B195D in detail 
based on their JfST/WFPC2 imaging survey of young 
massive GCs in M31. They initially selected VDBO- 
B195D as two YMCs in M31, but their WFPC2 images 
showed unequivocally that these two sample objects are, 
in fact, the same cluster. In addition, the HST images 
clearly confirmed that VDB0-B195D is a real cluster. 
However, it is difficult to establish whether it is more 
similar to ordinary open clusters, similar to those in the 
disk of the Milky Way, than to YM Cs that may evolv e 
to become disk GCs (see for details Perina et al.l l2009^. 
Sp ectral observations of VDB0-B195D were obtained 
by Ivan den Bergh I ()1969[ ) — yielding classi fication spec- 
tra an d the object's radial velocity — and iPerrett et al.l 
(|2002l) . who used them for determination of its radial 
velocity and metallicity. 

2.2. Archival images of the BATC Multicolor Sky 
Survey 

Observations of the YMC VDB0-B195D were obtained 
with the BATC 60/90cm Schmidt telescope located at 
the XingLong station of the National Astronomical Ob- 
servatory of China (NAOC). This telescope is equipped 
with 15 intermediate-band filters covering the optical 
wavelength range from 3000 to 10,000 A. The filter sys- 
tem was specifically designed to avoid contamination by 
the brightest and most variable night-sky emission lines. 
Descript ions of th e BATC photometric system can be 
found in lFan etalj (|1996f ). Before February 2006, a Ford 
Aerospace 2kx2k thick CCD camera was installed, with 
a pixel size of 15 /im and a field of view of 58' x 58', yield- 
ing a resolution of 1.7" pixel"^. Since February 2006, a 
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new E2V 4kx4k thinned CCD with a pixel size of 12 
/xm has been in operation, featuring a resolution of 1.3" 
pixel" ^. The blue quantum efficiency of the new, thinned 
CCD is 92.2% at 4000 A, which is much higher than for 
the old, thick device (see for details iFan et al.ll2009D . A 
field including VDB0-B195D in the a-c filters was ob- 
served with the thinned CCD, and in d-p bands with the 
thick CCD. Fig. 1 shows a finding chart of VDB0-B195D 
in the BATC g band (centered at 5795 A), obtained with 
the NAOC 60/9Gcm Schmidt telescope. We adopt an 
aperture with a radius of 15" (shown in Fig. 1) for the 
integrated photometry discussed in this paper. 

The BATC survey team obtained 61 images of VDBG- 
B195D in 15 BA TC filters betwee n January 2004 and 
November 2006. iFan et al.l ()2OO90 performed the data 
reduction of these images, which formed part of their 
M31-7 field. Table 1 contains the observation log, includ- 
ing the BATC filter names, the central wavelength and 
bandwidth of each filter, the number of images observed 
through each filter, and the total observing time per ffi- 
ter. Multiple images through the same filter were com- 
bined to improve image quality (i.e., increase the signal- 
to-noise ratio and remove spurious signal). 

2.3. Intermediate-band photometry of VDB0-B195D 

We determined the intermediate-band magnitudes of 
VDB0-B195D on the combined images using a standard 
aperture photomet ry app roach, i.e., the phot routine in 
DAOPHOT (Stct sonill987t l. Calibration of the magnitude 
zero level in the BATC photometric system is similar 
to that of the spectrophotometric AB m agnitude sys- 
tem. For flux calibration, the Oke-Gunn ()Oke fc GunnI 
HgSS") primary flux standard stars HD 19445, HD 84937, 
BD -|-26°2606, and BP -H7°4708 we re observed during 
photometric nights (|Yan et al.l [20001) . VD B0-B195D is 
located in the M31-7 field of iFan et all (|2009[) . The 
absolute flux of the M31-7 field was calibrated based 
on secondary standard transformations using the M31- 
1 field, which was calibrated, in tur n, by the four Oke - 
Gunn primary flux standard stars bv lJiang" "ei~aI1([2003h . 
Since VDB0-B195D is an extended object, an appro- 
priate aperture size must be adopted to determine its 
total luminosity. The (radial) photometric asymptotic 
growth curve, in all BATC bands, flattens out at a ra- 
dius of ~ 15". Inspection ensured that this aperture 
is adequate for photometry, i.e., VDB0-B195D does not 
show any obvious signal beyond this radius. In addi- 
tio n, this aperture is nearly the same as that adopted 
by iPerina et al.l ()2009D to determine the cluster's pho- 
tometry in the BVRI bands, based on the M31 imaging 
survey of Massey (2006) (see §2.4 below). Therefore, we 
use an aperture with r sa 15" for integrated photome- 
try, i.e., r = 9 pixels for the 2kx2k thick CCD camera, 
and r = 12 pixels for the 4kx4k thinned CCD camera. 
VDB0-B195D is projected onto the disk of M31, where 
the background is bright and fluctuates, potentially as 
a function of distance from the cluster center. To avoid 
contamination from background fluctuations, we adopted 
annuli for background subtraction spanning between 10 
and 15 pixels for the 2kx2k thick CCD camera, and from 
13 to 20 pixels for the 4x4k thinned CCD camera, both 
corresponding to ~ 17-26". While these annuli are spa- 
tially as close as possible to the region dominated by 



cluster light (so that any differences in background flux 
are minimized), they are wide enough to average out any 
expected background fluctuations. The calibrated pho- 
tometry of VDB0-B195D in 15 filters is summarized in 
column (6) of Table 1, in conjunction with the la magni- 
tude uncertainties, which include uncertainties from the 
calibration erro rs of both the M 3 1-1 field standard stars 
(see for details iFan et al.l 120091: iJiang et al.l I2003D and 
'the secondary standard stars' in common between the 
M31-1 and M31-7 fields used for calculation of the mean 
magnitude offsets between the standard and in strumen- 
tal magnitudes (see for details iFan et all l2009| ) , as well 
as those resulting from our daophot application. 

2.4. Optical broad-hand and near-infrared 2MASS 
photometry of VDB0-B195D 

Four indepen dent sets of ph o tomet ric data exist for 
VDB0-B195D. Ivan den Ber"iiri ()1969[ ) obtained UBV 
photo metry using observation s of the 200-inch Hale tele- 
scope, iBattistini et al.l ()1987D performed UBVR pho- 
tometry based on photographic plates observed with the 
152 cm Ritchey-Chretien //8 telesc ope of t he Un iver- 
sity of Bologna in Loiano, King & LuptonI (|1991[ ) ob- 
tained UBV photometry for VDB0-B195D using obser- 
vations with the University of Hawaii's 2.2 m telescope 
on Mauna Kea using the //lO seco ndary and coronene - 
coated 584 x 416 GEC CCD, and ISharov et al.) ([I99l 
performed UBV photometry based on photo-electric ob- 
servations with the 2.6 m Shain telescope of the Crimean 
Astrophysical Observatory. In addition, in the Revised 
Bologna Cat a logue (RBC) of M31 GCs published by 
iGalleti et ail (|2004[ ). the photometric data of VDBO- 
B195D in optical bands are b ased on IBattistini et al.l 
(I1987D ISharov e t al." (1995") , and transformed to the 
reference system of [Barmby et al. (200_Q) by applying off- 
sets derived from objects in com mon between the rele - 
vant catalog and the data set of iBarmbv et al.l ()2000D . 
In the RBC, VDB0-B195D was regarded as two objects. 
We list these photometric data in Table 2 for comparison. 
Note that, in the latest RBC incarnation (version 3.5, up- 
dated on 27 March 2008), VDB0-B195D is included as a 
si ngle object. 

IGalleti et al] (|20Q1 also determined 2MASS JHK^ 
photometric magnitudes for VD B0-B195D (tr ansformed 
to the CIT photometric system: IsFas et al.| [T982. 1983), 
which we have in cluded in Table 3. In addition, 
IPerina et aH (|2009( ) realized that VDB0-B195D is a par- 
ticularly extended object and t hat it is possible that 
the photometry of ISharov et al.l (|1995f) (compiled in the 
RBC) was obtained with apertures that were not large 
enough to include all of its flux. Therefore, they redeter- 
mined its photometric values i n the BV RI ba nds based 
on the M31 imaging survey of lMassevI jlOOe') using an 
aperture with r = 14.4", which are also listed in Table 3. 

From a comparison of the valu es in Tables 2 and 3, it 
is clear that the magnitudes of Ivan den Bergh I (Il969l) 
are brighter, while the results of the three other ref- 
erences are consist ent. The magnitudes determined by 
IPerina et all (|200l ) are much brighter, however, because 
of their careful inclusion of all of the cluster's flux. To 
compare our photometric results with previously pub- 
lished values, we transformed the magnitudes of VDBO- 
B195D in the BATC intermediate bands to broad-band 
U i?V^i?/-equivalent photometry based on the relation- 
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Fig. 1.— Image of VDB0-B195D in the BATC g band, obtained with the NAOC 60/90cm Schmidt telescope. VDB0-B195D is circled 
using an aperture with a radius of 15". The field of view of the image is 11' X 11'. 



ships obtained by IZhou et"aLl ()2003D . These are also 
Hsted in Table 3, and the uncertainties include those 
originating fr om the tra nsformation based on the rela- 
tionships of Z hou et al.l (2003) and their calibration er- 
rors (column 5 of their Table 3). In Fig. 2, we show the 
result of the comparison. In general, the oth er photomet- 
ric da ta are fainter than ours and those of iPerina et al.l 
Hooi). Fig. 2 and Table 3 show that our new VRI mag- 
nitudes agree with the results of Perina et al.l |2009), and 
that the B magnitude obtained in this pap er is 0.32 mag 
brighter than that of IPerina et alj (I2009D. Consid ering 
the photometric errors of both IPerina et all (|2009f ) and 
our current study, these two -B-band photometric results 
are consistent within 2a. In addition, we should keep 
in mind that, although the VRI magnitude s obtained in 
this p aper are consistent with the results of IPerina et al.l 
(|2009f) within la, the disagreement in B magnitudes at 
this level is understandable. This is caused by the fact 
that the original photometry in the present paper was ob- 



tained in the proprietary BATC filters and transformed 
to the UBVRI sy stem using transformation equations. 
I Zhou et"al] (|2003l ) determined these conversions based 
on the b r oad-band U BVR I magnitudes of 48 sta r s from 
iLandoltj (fl98l [1991 and iGaladi-Enn'guez et al.l pOOOl ) 
in the Landolt SA95 field, and their photometric data 
in the 15 BATC intermediate-band filters. In addition, 
the central wavelengths and bandwidths of the BATC 
and U BVRI systems differ. In fact, a similar significant 
disagreement of i3-ba nd photomet r ic dat a for some M31 
GCs was reported bv iWang et al.l (|2010f ). citing similar 
arguments. 

3. STELLAR POPULATION OF VDB0-B195D 

3.1. Stellar populations and synthetic photometry 

To determine the age and mass of VDB0-B195D, 
we compared its SED with theoretical stellar popula- 
tion synthesis models. The SED consists of photo- 
metric data in the 15 BATC intermediate bands ob- 
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R (BATC) / (BATC) 



Fig. 2. — Comparison of photome tric data from diffe rent sources with new determinations in this paper for VDB0-B195D. The data 
points shown as black dots are from Peri na et ahl l|2009l ) . 



tained in this paper and optical broad- band BVRI 
and 2 MASS near- infrared JHKs data from lPerina et all 
(I2009D. listed in Table 3. We used the galev 



SSP models (e.s.. IKurth et al.lll999l: 


Schulz et alJ 120021: 


lAnders & Fritze-v. AlvenslebenI 12003 


) for our compar- 



isons. The GALEV SSPs are based on the Padova stellar 
isochrone s, with the most recent versions using the up- 
dated iBertellOr^ (|1994| ) isochrones (which include the 
thermally pulsi ng asymptotic giant-branch phase), and a 
ISalpeteii (jl955D stellar IMF with lower- and upper-mass 
limits of 0.10 and between 50 and 70 Mq, respectively, 
depending on metallicity. The full set of models spans 
the wavelength range from 9lA to 160 fim. These mod- 
els cover ages from 4 x 10^ to 1.6 x 10^" yr, with an age 
resolution of 4 Myr for ages up to 2.35 Gyr, and 20 Myr 
for greater ages. The galev SSP models include five 
initial metallicities, Z = 0.0004,0.004,0.008,0.02 (solar 
metallicity), and 0.05. 

Since our observational data consist of integrated lumi- 
nosities through the set of BATC filters, we convolved the 
GALEV SSP SEDs with the BATC intermediate-, opti- 
cal broad-band BVRI, and 2MASS filter-response curves 
to obtain synthetic optical and near-infrared photometry 
for comparison. The synthetic i"^ filter magnitude can 
be computed as 



-2.5 log -— 48.60, 



(1) 



where F,y is the theoretical SED and (pi the response 
curve of the i"^ filter of the BATC, BVRI, and 2MASS 
photometric systems. Here, varies with age and 



metallicity. Since the observed magnitudes in the BVRI 
and 2MASS photometric systems are given in the Vega 
system, we transformed them to the AB system for our 
fits. 

3.2. Reddening and metallicity of VDB0-B195D 

To obtain the intrinsic SED of VDB0-B1 95D, its pho- 
tomet ry must be dereddened. To date, onlv lPerina et al.l 
()2009f ) obtained reddening values for VDB0-B195D. 
They compared the observed CMD with theoretical 
isochrones and determ ined E{B — V) — 0.20 ± 0.03 mag. 
ICaldweU et all (l2009h were unable to derive the cluster's 
reddening value because of the presence of a foreground 
field star, so they adopted E{B -V) ^ 0.28 ± 0.17 mag 
(external rms error), equivalent to the mean reddening 
of the young clusters in M31. In this pa per, we ther efore 
adopt the reddening value from iPerina et al.l (j2009D . 

In addition, cluster SEDs are afii'ected by age and 
metallicity effects. Therefore, we can only accurately 
constrain a c luster 's age if the metallicity is known. 
IPerina et all ()2009[ ) found that the CMD of VDBO- 
B195D, based on their ffS'T/WFPC2 observations, is 
best reprodu c ed b y the solar-metallicity models of 
iGirardi et al.l (|2002| ). We therefore adopt solar metal- 
hcity for VDB0-B195D. 



The 



3.3. The 'lowest-luminosity-limit' test 

lowest-lumino sity limit (LLL; 



iCervino fc Luridianal I2004D implies that it is mean 
ingless to compare a cluster with population synthesis 
models to obtain its age and mass if its integrated 
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luminosity is lower than the luminosity of the most 
luminous star included in the model for the relevant 
age. The LLL method states that clusters fainter than 
this limit cannot be analyzed using standard procedures 
such as minimization of the observed values with 
respe ct to the mean SSP models (see also I Barker et alJ 
|2008() . Below the LLL, cluster ages and masses cannot 
be obtained self-consistently. To take into account the 
effects on the integrated luminosities of statis tically 
samp l ing the stellar I MF (e.g., iCervifio et aH 120001 
ervino fc Luridi ana 2004!), we used the theoret- 
ical Padova isochrones at http://stev.oapd.inaf.it/cmd 
(CMD2.2). This interactive Web interface provides 
isochrones for a number of photometric systems, in- 
cluding optical broad-band, 2MASS, and the BATC 
data used here. We obtai ned the sola r- metal licity 
(Z = 0.019) isochrones of iMarigo et aT (200 81), as 
recommended by CMD2.2, based on the (iSalpeten ll955| ) 
IMF so as to match the IMF selection for our age and 
mass determinations of VDB0-B195D in §3.4 (see §3.1 
for details). 

Figure 3 shows the LLL values as a function of age for 
the different filters used in this paper. These luminosi- 
ties are obtained by identifying the most luminous star 
on each isochrone for the relevant passband. The gray 
area shows the cluster's absolute luminosity, assuming a 
distance modulus of (m — M)n = 24.47 mag (785 kpc) for 
M31 (Mc Connachie et all [20051 ). The upper luminosity 
limit has been corrected for extinction, based on a red- 
dening value of E{B — V) = 0.20 mag. The interstella r 
extinction curve. Ax, is taken from lCardelli et all (|1989[ ). 
Ry = Av/E{B -V)= 3.1. 

We see that VDB0-B195D does not lie below the LLL 
in any of the passbands used here. This means that, in 
general, VDB0-B195D can host the most luminous star 
that would be present theoretically for the given age of 
the cluster. 

3.4. Fit results 

In the previous section, the LLL test proves that the 
luminosity of VDB0-195D is higher than the luminosity 
of its brightest star expected for a given cluster age, i.e., 
that using SSP models is not completely meaningless. In 
addition, the bright absolute magnitude of VDB0-195D 
allows us to consider a possibility that the cluster is mas- 
sive enough and IMF sampling effects should not strongly 
impact the fitting results. So we will determine the clus- 
ter's age and mass estimates based on direct compar- 
isons with SSP mean values in this section. However, we 
should keep in mind that this approach is a compromise. 
In fact, the fitting results (Fig. 4 and Table 5) show 
probable problem even for relative massive clusters. 

We use a minimization test to determine which 
GALEV SSP models are most compatible with the ob- 
served SEDs, 

^ =2. -2 ' (2) 

i=l » 

where m™°'^(i) is the integrated magnitude in the i"' 
filter of a theoretical SSP at age t (for solar metallicity), 
m™**^ is the intrinsic, integrated magnitude, and ai is the 



magnitude uncertainty, defined as 

'^i = ^lhs,i + ^mod.i + * CTrcd)^ + CTmd,i- (3) 

Here, (Tobs.i is the observational uncertainty from col- 
umn (6) of Table 1 and column (2) of Table 3, (Tmod.i 
is the uncertainty associated with the model itself, di-cd 
is the uncertainty in the reddening valu e, and R\, — 
A\JE (B — V), where A\. is taken from iCardelli et aTl 
(|1989D . Rv = Av/E{B - V) = 3.1, and cr^d.^ is the 
uncertainty in the d istance modulus, for the i*^ filter. 
I Chariot et al.l (1996) estimated the uncertainty associ- 
ated with the term dmod.i by comparing the colors ob- 
tained from different ste llar evolut i onary track s and spec- 
tral libraries. Following fMa et al.l (|2007l l2009( ) . we adopt 
= 0-05 mag. 

iPerina et al.] p009l ) pointed out that VDB0-B195D is 
a particularly extended object and that the photometric 
measurements of van den Bcrgh (1969), Battistini et al] 
(|1987f) . IKing~FLupton (1991), and .Sharov et al.. (,19951 ) 
did not include all of its flu x The refore, we adopt the 
photometry of IPerina et al.l ()2009D to fit the observed 
SED with theoretical SSPs for our age determination. 
The fit yielding the minimum value (x^(min)) was 
adopted as the best fit and we adopted the correspond- 
ing age value, 60.0 ± 8.0 Myr. In addition, our best- 
fitting age estimate of 60.0 ± 10.0 Myr results from us- 
ing the (redder) k-p and IJHKs photometry; using only 
the blue part of the cluster's SED [B, a-e, where any 
effects caused by stochasticity may be smaller) yields 
an age of 72.0 ± 34.0 Myr. The uncertainty was esti- 
mated using confidence limits. If < X^i^^^)/^ + 
the resulting age is within the 68.3% probability range; 
here, = 21 is the number of free parameters, i.e., the 
number of observational data points minus the number 
of parameters used in the theoretical model. Therefore, 
the accepted age range is derived from those fits that 
have x^(min)/j^ < X^/'^ < X^(inin)/i/ + 1. The best 
reduced-x^ — defined as x^(niiii) = x^(nrin)/i/ — and age 
are listed in Table 4. The best fit to the SED of VDBO- 
B195D is shown in Fig. 4, where we display the intrinsic 
cluster SED (symbols with error bars), as well as the 
integrated SED (open circles) and spectrum of the best- 
fitting model. From Fig. 4, we note that the observa- 
tional data in the b, d, o, and p BATC filters and in the 
Ks band do not match the best-fitting model very well 
(the difference is approximately 0.3 mag). Photometric 
uncertainties in these filters may cause some differences, 
although this might not be the main reason for the dis- 
crepancy. As we know, observational star clusters' SEDs 
are affected by age, metallicity and reddening. If the 
reddening value and metallicity adopted in this paper 
are not problematic, discrepancy between our observa- 
tions and the best-fitting model may refiect the difficulty 
in achieving an appropriate (but formal) fit of an SED of 
a single, real cluster by SSP models. However, as we will 
see below, the reddening value adopted in this paper may 
be bigger than the actual reddening of VDB0-B195D. In 
addition, the differences between the photometric data 
and the model in Fig. 4 show a somewhat systematic 
behavior with wavelength: in bluer passbands the clus- 
ter seems to be more luminous than predicted by the 
model, while in redder passbands it is fainter than the 
corresponding model predictions. A blue excess and red 
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Fig. 3. — Lowest-luminosity limit for the filters used in this paper. The curves indicate the luminosities of the most luminous star on each 
isochrone for the relevant passb and . The light-gray area shows the absolute magnitudes of VDB0-B195D based on a reddening value of 
E(B - V) = 0.20 mag (Perina ct 'alll2009l ) . We used a distance modulus of (m - M)o = 24.47 mag (785 kpc) for M31 (,McConnachie et al.l 
120051) to calculate the absolute magnitudes. 
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deficiency in the observed SED with respect to the model 
predictions may indicate a shortage of red giants (RGs), 
which can occur when the cluster is either younger or 
less massive (or both) than the corresponding best-fitting 
model suggests. In other words, IMF discreteness may 
play a role: due to a relatively longer main-sequence 
(MS) phase and shorter RG phase, a random young clus- 
ter is typically bluer than predicted by SSP models. At 
the same time, we find that the reddening value adopted 
afl^ects the fitting result greatly. In fact, the best fit to the 
SED of VDB0-B195D improves a great deal when adopt- 
ing a smaller reddening value such as E{B ~ V) = 0.1: 
X^(min) = 0.73; the resulting age (64.0 ± 8.0 Myr) is 
nearly the same as one (60.0 ± 8.0 Myr) obtained with 
E{B -V)= 0.2. 

We next determined the mass of VDB0-B195D. The 
GALEV models include absolute magnitudes (in the Vega 
system) in 77 fihers for SSPs of 10^ Mq, includi ng 66 fil- 
ters o f the HST, Johnson U BVRI (s ee for detailslLandoltl 
Il983f) . Cousins RI (see for details iLandoltl |1983D, and 
JHK (jBesseh fc Br ett 1988) systems. The difference be- 
tween the intrinsic absolute magnitudes and those given 
by the model provides a direct measurement of the clus- 
ter mass, in units of 10^ Mq. However, we should keep 
in mind that this is only correct for cluster masses above 
10^ Mq. We estimated the mass of VDB0-B195D using 
magnitudes in all of the BVRI and JHKs bands. There- 
fore, we transformed the 2MASS JHKs magnitudes to 
the photometric syste m ofjB esscll & Brett (1988) using 
the equations given bv ICarpenter. (2001). The resulting 
mass determinations for VDB0-B195D are listed in Ta- 
ble 5 with their \a uncertainties including contributions 
from uncertainties in extinction and distance modulus. 
From Table 5, we see that the mass of VDB0-B195D ob- 
tained based on the magnitudes in different filters is very 
different. (The highest mass obtained, based on the B- 
band magnitude, is 0.5 x 10^ Mq more massive than 
that obtained using the Kg magnitude.) In addition, the 
mass estimates differ systematically with filters. Pro- 
vided that VDB0-B195D is massive enough to be fitted 
by SSP models, a systematic trend of masses based on 
different passbands may indicate a problem with redden- 
ing value adopted for the cluster. If the actual reddening 
is smaller than the adopted value, the actual luminosity 
would be overestimated. This effect is small in redder 
filters but strong in bluer filters. As discussed in age es- 
timation, a smaller reddening value can improve the fit- 
ting result greatly. In fact, a smaller reddening value can 
reduce the mass discrepancies based on the magnitudes 
in different filters. When we adopted E{B — V) = 0.1, 
the mass of VDB0-B195D obtained based on the magni- 
tudes in different filters is the same within Icr. We list 
these estimates in Table 7. From Table 5, we know that 
the mass of VDB0-B195D obtained in paper is between 
(1.1 — 1.6) X 10^ Mq when the reddening value is adopted 
to be {B-V)^ 0.2. 

4. SUMMARY AND DISCUSSION 

VDB0-B195D was previously shown to be a massive 
cluster based on ffS'T/WFPC2 observations. Its color is 
extremely blue and it is very bright, particularly in blue 
bands. In addition, VDB0-B195D is an extended object, 
and most previous photometric measurements did not 
include its full flux distribution (see Perina ct al.. .20091 . 



for details). 

In this paper, we obtained the cluster's SED in the 
15 BATC intermediate-band filters. We subsequently 
determined its age and mass by comparing our multi- 
color photometry with theoretical stellar population syn- 
thesis models. Our multicolor photometric data con- 
sist of 15 intermediate-band filters obtained in this pa- 
per, and bro ad-band BVRI and 2MASS JHK^ from 
iPerina et al.) (|2009), covering a wavelength range from 
3000 to 20,000 A. Our rcsuffs show that VDB0-B195D is 
a genuine YMC in M31. 

To under st and the real nature of the BLCCs, 

IPerina et al.) pOOl [2010l ) performed an HST imag- 
ing s urvey of 20 BLCC s in MSl's disk. As a test 
case, IPerina et al.l (j2009D presented details of the data- 
reduction pipeline that will be applied to all survey data 
and describe its application to VDB0-B195D. They es- 
timated the object's age, by comparison o f the observed 
CMD with theoretical isochrones from iGirardi et al.l 
(|2002f) . at ~ 25 Myr. In addition, they constrained re- 
alistic upper and lower limits to the cluster's age, in- 
dependent of the adopted metallicity, within the rela- 
tively narrow range from 12 to 63 Myr. Usiiig Mar aston's 
SSP models of solar metallicit y (lMarastonl ll99l [2005), 
iSalpeteii (|1955D and iKroupal daOOl") IMFs, and photo- 
metric values in the V and 2MASS J, H, and A's bands, 
IPerina et ahl (|2009D concluded that the mass of VDBO- 
B195D is > 2.4 x 10'* Mq, with their best estimates in 
t he range - (4 - 9) x l O"^ Mq. 

ICaldwell et al.l (|2009D presented an updated catalog of 
1300 objects in M31, including spectroscopic and imag- 
ing surveys, based on images from either the LGGS or 
the DSS and spectra taken with the Hectospec fiber po- 
sitioner and spectrograph on the 6.5 m MMT. They de- 
rived ages and reddening values for 140 young clusters 
by comparing their observed spec tra with model spectr a 
from the Starburst99 SSP suite (jLeitherer et al.lll999D . 
The results show that these clusters are less than 2 Gyr 
old, while most have ages between 10^ and 10^ yr (the age 
of VD B0-B195D they derive is logage/yr — 7.6). In ad- 
dition, ICaldwell et al.l ()2009[ ) also estimated the masses 
of these young clusters usi ng F-band photomet ry and 
model mass-to-light ratios (jLeitherer et al.|[l999( ) corre- 
sponding to the derived spectroscopic ages. This resulted 
in masses ranging from 2.5 x 10^ to 1.5 x 10^ Mq. The 
mass of VDB0-B195D obtained bv ICaldwell et all (|2009[ ) 
is logAfciA^o =5.1 (no uncertainty quoted). 

We compare the various age and mass estimates of 
VDB0-B195D in Table 6. Ou r newly obta i ned a ge is 
older than the e stima tes of both IPerina et al.l (j2009() and 
ICaldwell et all (|2009t) . while the mass obtained in this 
paper is higher than the estimate of IPerina et al.l (I2009D 
and co nsistent with the determination of ICaldwell et al.l 
(|2009| ) . However, our r esults are ii i agreement wit h those 
of both IPerina et al.) (|2009l ) and ICaldwell et all (|2009l ) 
within 3(7. The age and mass obtained in this paper 
confirms that VDB0-B195D is genuinely a YMC in M31. 

As we know, SSP models describe a very special case of 
a continuous distribution of stellar mass (or light) along 
isochrones. This is well approximated by clusters with 
masses larger than 10^ Mq. Also, for cluster masses of 
about 10^ Mq, SSP models can probably still be applied 
since a systematic difference between SSP models and ob- 
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Fig. 4. — Best-fitting, integrated tlieoretical GALEV SEDs compared to the intrinsic SED of VDB0-B195D. Tfie pfiotometric measurements 
are shown as symbols with error bars (vertical for uncertainties and horizontal for the approximate wavelength coverage of each filter). 
Open circles repres ent the calculated magnit udes of the model SED for each filter. We used a distance modulus of (m — M)o = 24.47 mag 
(785 kpc) for M31 IIMcConnachie et al.ll2005l ) to calculate the absolute magnitudes. 



servations should, on average, be smaller than 0.05 mag 
for clusters older than 10 Myr (see Fig. 3 in Piskunov 
et al. 2009). However, from the results of this paper, 
we may conclude that, probably, a formal fitting of SSP 
models to observed SEDs cannot be used without cau- 
tion even for relatively massive (or apparently massive) 
clusters, and it is highly doubtful that this approach can 
be applied in a routine work providing accurate cluster 
parameters. The relative accuracy of 10% for age and 
20% found for the mass of VDB0-B195D seems to be 
rather formal and not very confident. In addition, obser- 
vational star clusters' SEDs are affected by reddening, an 
effect that is also difficult to separate from the combined 
effects of age and metallicity (Calzetti 1997; Vazdekis et 
al. 1997; Origlia et al. 1999). Only the metalhcity and 
reddening are derived accurately (and, ideally, indepen- 
dently), these degeneracies are largely (if not entirely) 
reduced, and ages can then also be estimated accurately 
based on a comparison of multicolor photometry span- 
ning a significant wavelength range (de Grijs et al. 2003b; 



Anders et al. 2004) with theoretical stellar population 
synthesis models. It is true that the discrepancy between 
our observations and the best- fitting model is great, and 
the mass of VDB0-B195D obtained based on the magni- 
tudes in different filters is very different. However, when 
we adopt a smaller reddening value, the results improve 
greatly. So, we conclude that the actual reddening value 
of VDB0-B195D may be smaller than E{B -V)= 0.2. 
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TABLE 1 

BATC PHOTOMETRY OF THE M31 YMC VDB0-B195D. 



i-'ilter 


Central wavelength 
(A) 


Bandwidth 
(A) 


Number oi images 


Exposure time 
(hours) 


Magnitude 


a 


3360 


222 


6 


2.0 


15.51 ± 0.14 


b 


3890 


187 


6 


2.0 


14.73 ± 0.11 


c 


4210 


185 


3 


0.8 


14.70 ±0.07 


d 


4550 


222 


3 


1.0 


14.49 ± 0.10 


e 


4920 


225 


3 


1.0 


14.67 ±0.05 


f 


5270 


211 


3 


1.0 


14.65 ± 0.05 


9 


5795 


176 


3 


1.0 


14.59 ± 0.04 


h 


6075 


190 


3 


1.0 


14.56 ± 0.02 


i 


6660 


312 


3 


1.0 


14.50 ±0.02 


3 


7050 


121 


5 


1.7 


14.51 ± 0.06 


k 


7490 


125 


3 


1.0 


14.47 ±0.05 


m 


8020 


179 


3 


1.0 


14.43 ± 0.07 


n 


8480 


152 


6 


2.0 


14.49 ± 0.05 





9190 


194 


6 


2.0 


14.49 ± 0.05 


P 


9745 


188 


6 


2.0 


14.50 ±0.05 



TABLE 2 

Comparison of broad-band photometry of VDB0-B195D. 

Filter Mag" Mag" Mag'' Mag" Mag° MagJ 

U 14.66 15.11 15.12 ±0.012 14.97 ± 0.01 15.110 15.140 

B 15.14 15.39 15.49 ±0.010 15.31 ± 0.01 15.410 15.510 

V 14.94 15.19 15.32 ±0.013 15.06 ± 0.01 15.190 15.280 

_R, 15.27 14.920 



Ivan den Berghl 119691 ): ' iBattistini et all 119871) : IKing fc LuptonI 119911) . uncertainties are the median 
uncertainties in the mean for all sample cluster measurements; ISharov et al.l 119951 ): Photometry from [Ualleti et all 120041) , based on 
IBattistini eFaLl | |T987|) : •''Photometry from lGalleti et al.l ll200l ), based on ISharov eFaLl | |T995|) . 



TABLE 3 

Recently determined photometry' for VDB0-B195D. 



Filter 


Mag" 


Mag" 


Mag'' 


U 






14.37 ±0.22 


B 


14.94 ±0.09 




14.62 ±0.13 


V 


14.67 ±0.05 




14.67 ±0.05 


R 


14.45 ±0.11 




14.60 ±0.06 


I 


14.01 ±0.11 




14.19 ±0.10 


J 


13.26 ±0.07 


13.78 ± 0.03 




H 


12.76 ±0.12 


13.15 ± 0.04 




Ks 


12.77 ±0.15 


12.96 ± 0.03 





" IPerina et aP ifSOOgP : tealleti et all l(200l 1: ''This paper. 
TABLE 4 

Age estimate of VDB0-B195D based on the the galev models. 



Age 


log (Age) 


X^(min) 


(Myr) 


[yr] 


(per degree of freedom) 


60.0 ±8.0 


7.78 ±0.05 


2.2 



TABLE 5 

Mass estimates (and uncertainties) of VDB0-B195D based on the galev models. 



B 


V 


R 


I 


J 


H 




Mass (10^ Mq) 


1.6 ±0.18 


1.6 ±0.13 


1.4±0.17 


1.4 ±0.16 


1.4 ±0.13 


1.3 ±0.16 


1.1 ±0.17 



Ma et al. 



TABLE 6 

Comparison of age and mass estimates of VDB0-B195D. 



Age^ 


Age^ 


log (Age)=' 


log (Age) 2 


Mass^ 


Mass^ 


log (Mass)^ 


log (Mass)^ 


(Myr) 


(Myr) 


[yr] 


[yr] 


(10* Mo) 


(105 Mq) 




[Mq] 


25 


60.0 ±8.0 


7.6 


7.78 ± 0.05 


4-9 


1.1 - 1.6 


5.1 


5.0 - 5.2 



J |Perina et al] ||2009|) : ^This paper; ' ICaldwell et aO l|2009l) . 



TABLE 7 

Mass estimates (and uncertainties) of VDB0-B195D based on the galev models with E(B — V) = Q. 



B 


V 


R 


I 


J 


H 


Ks 


Mass {Vf> Mq) 


1.1 ±0.12 


1.2 ±0.10 


1.1 ±0.14 


1.2 ±0.14 


1.3 ±0.12 


1.2 ±0.16 


1.1 ±0.17 



